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ABSTRACT 

We present black hole masses and accretion rates for 182 Type 1 AGN in COSMOS. We estimate 
masses using the scaling relations for the broad H/3, Mgn, and Civ emission lines in the redshift 
ranges 0.16 < z < 0.88, 1 < z < 2.4, and 2.7 < z < 4.9. We estimate the accretion rate using an 
Eddington ratio Lj /L Edd estimated from optical and X-ray data. We find that very few Type 1 AGN 
accrete below Lj/LEdd ~ 0.01, despite simulations of synthetic spectra which show that the survey 
is sensitive to such Type 1 AGN. At lower accretion rates the BLR may become obscured, diluted 
or nonexistent. We find evidence that Type 1 AGN at higher accretion rates have higher optical 
luminosities, as more of their emission comes from the cool (optical) accretion disk with respect to 
shorter wavelengths. We measure a larger range in accretion rate than previous works, suggesting 
that COSMOS is more efficient at finding low accretion rate Type 1 AGN. However the measured 
range in accretion rate is still comparable to the intrinsic scatter from the scaling relations, suggesting 
that Type 1 AGN accrete at a narrow range of Eddington ratio, with Lj/LEdd ~ 0.1. 
Subject headings: galaxies: active — quasars 



1. INTRODUCTION 

Supermassive black h oles (SMBHs) reside i n al- 
most all local galaxi es (jKormendv fc Richstond 119951 : 
iRichstone et all Il998l ). The mass of the SMBH is 
observed to be tightly correlated with the mass, lu- 
minosity, and velocity disp ersio n of the host galaxy 
bulge (e.g..lMagorrian et al.l ll998t lGebhardt etallboori : 



iFerrarese fc MerrittJ I2000D . SMBHs grow by accre- 
tion as active galactic nuclei (AGN), and all mas- 
sive gal a xies have one or more o f the se active phases 
HSolt'arj 119821 : iMagorrian et alj 119981 : IMarconi et all 
|2004[ ). More lu minous AGN are observed to peak 
at higher redshift (lUeda et al.ll2003t iBrandt fc Hasingerl 
120051 IBongior no et al.|[2007| ) . exhibiting "downsizing" by 
analogy to the preference of luminous and massive galax- 
ies to form at high redshifts. Both downsizing and the 
correlations between SMBH and the host bulge suggest 
that the growth of AGN and the forma tion of galaxies are 
directly connected through feedback (|Silk fc Reeslll998l : 
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iDi Matteo. Springel. fc Hernquistlf2005f ). 

Understanding the role of AGN in galaxy evolution re- 
quires measurements of SMBH mass and accretion over 
the cosmic time. The SMBH mass can be directly esti- 
mated by modeling the dynamics of nearby gas or stars, 
but this requires high spatial resolution and is limited 
to HST observations of nearby galaxies. Reverberation 
mapping of Type 1 AGN (with broad emission lines) 
uses the time lag between variability in the continuum 
and the broad line region (BLR) to estimate the radius 
of the broad line region, R blb. = ctiag (for a review, 
see iPeterson fc Ben tz 2006) . Then, if the broad line re- 
gion virially orbits the source of the continuum emis- 
sion, the SMBH mass is Mbh = J ' RBLKV 2 f whml where 
/ represents the unknown BLR geometry and Vf w h m is 
the velocity width of the broad emission lin e. This tech- 
nique has many pote ntial systematic errors (|Krolikll2001l : 
IMarconi et al.ll2008f) . but its mass e stimates agree with 
those from dynam ical estimators (jDavies et al.l 120061 : 
lOnken et al.ll2007l) and those from the Mbh -v* corre- 
lation (|Onken et al.|[200l iGreene fc Hoi [20061. In prin- 
ciple, reverberation mapping can be applied to AGN at 
any redshift, but in practice, the need for many periodic 
observations has limited reverberation mapping mass es- 
timates to only ^35 local AGN. 

Instead, the vast majority of AGN mass estimates 
have come from a set of scaling relations. Reverbera- 
tion mapping data led to the discover y that Rblr corre - 
lates with the continuum luminosity dKaspi et al. 
with Rblr 1 ~ L ", where a ~ 0.5 (|Bentz et al. 
iKaspi et al.ll2007T ). This allows for estimates of Mbh 
from single epoch spectra with scaling relations: 
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Some authors replace the conti nuum luminosit y with 
the recombination line luminosity (iWu et al. 2004) or th e 
FWHM with the second moment a (j Collin et aiKl2006ft . 
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yielding minor sys tematic differences in estimated Mbh 
(|Shen et al.1 12008ft . While these scaling relations are 
based upon reverberation mapping of only local AGN, 
the method is based upon the abi lity of the central e ngine 
to ionize the broad line region (jKaspi et all 12000). and 
there is no physical reason to sugges t that the ionization 
of AGN should evo l ve wit h redshift (|Dietrich fc Hamannl 
120041: IVestergaardl |2004| ). Thus the scaling relations 
can be used to study t he distribution a n d evo lution of 
Type 1 AGN masses. iKollmeier et~aTI (|2006l) showed 
that Type 1 AGN tend to accrete at a narrow range 
of Eddington ratio, t ypically O.OlLEdd < L < Lsdd- 
IKollmeier etaLl ((2006) suggest a minimum accretion rate 
for Type 1 AGN, with AGN of lower accretion rate ob- 
serve d as "naked" Type 2 AGN without a broad lin e re- 
gion (|Hopkins et al.ll2008f) . iGavignaud et all (j2008ft ad- 
ditionally suggest that lower luminosity Type 1 AGN ac- 
crete less efficiently than brighter quasars. 

In this work we report black hole masses and study 
the demographics of 182 Ty pe 1 AGN in the C osmic 
Evolution Survey (COSMOS. IScoville et al.ll200l . We 
introduce the data and outline our spectral fitting in §2, 
and discuss the black hole masses, their associated errors, 
and our completeness in §3. We close with discussion of 
Type 1 AGN accretion rates in §4. All luminosities are 
calculated using h — 0.70, £1m = 0.3, f^A = 0.7. 

2. OBSERVATIONAL DATA 
2.1. Sample 

The Cosmic Evolution Survey (COSMOS, 
IScoville et all [20p) is a 2 deg 2 HST/ACS survey 
(|Koekemoer et all 120071 ) with ancillary deep multiwave- 
lcngth observations. The depth of COSMOS over such 
a large area is particularly suited to the study of low- 
density, rare targets like active galactic nuclei (AGN). 
The most efficient way to select AGN is by their X-ray 
emission, and XMM-N ewton observations of COSMOS 
(jCappelluti et alj I2008D reach fluxes of 1 x 10~ 15 cgs 
and 6 x 10" 15 cgs in the 0.5-2 keV and 2-10 keV bands, 
respectively. The matching of X-ray p oint sources 
to op tical counterparts is de scri bed by IBrusa et al.l 
(|2007l ) and IBrusa et alj (|2009l ). iTrump et alT ~(2009) 
performed a spectroscopic survey of XMM-selected AGN 
in COSMOS, revealing 288 Type 1 AGN with X-ray 
emission and broad emission lines in their spectra. Here 
we investigate a chief physical property, the black hole 
mass, for th e Type 1 AGN in th is survey. 

From the ITrump et al.1 (|2009l ) sample, we choose the 
182 Type 1 AGN with high-confidence redshifts and with 
H/3, Mgn, or Civ present in the observed wavelength 
range. That is, we select only z coa { > 3 AGN, empirically 
determined by ITrump et all (|2009f) to be at least 90% 
likely to have the correct classification and redshift. All 
of the Type 1 AGN spectra are dominated by blue power- 
law continua, with no obvious (beyond the noise) absorp- 
tion line signature from the host galaxy. The broad emis- 
sion line requirement restricts us to the redshift ranges 
0.16 < z < 0.88, 1 < z < 2.4, and 2.7 < z < 4.9. At 
these redshifts the AGN spect roscopy is > 90% c omplete 
to 4b < 22 ( see Figure 13 of lTrump et~aTll2009f) . where 
iX B is the AB magnitude from the COSMOS CFHT ob- 
servations. Our ability to measure a broad line width, 
however, is a slightly more complicated function of the 



spectral signal to noise. We characterize our complete- 
ness as a function of broad line width and S/N in §3.2. 

The majority (1 33) of the 182 AGN have spectra from 
Magellan/IMACS (jBigelow et al.lll998f) . with wavelength 
coverage from 5600-9200A and — 10A resolution. The re- 
maining bright 49 quasars have publicly available spec- 
tra from the Sloan Digital Sky Survey (SDSS) quasar 
catalog ([Schneider et al.ll2007l ). with 3800-9200A wave- 
length coverage and — 3A resolution. The resolution of 
both surveys is more than sufficient for measuring broad 
emission line widths; our narrowest broad emission line 
is 1300 km/s wide, compared to the resolution limits of 
IMACS (-600 km/s) and the SDSS (-200 km/s). 

2.2. Spectral Fitting 

The optical/UV spectrum of a Type 1 AGN can be 
roughly characterized as a power-law continuum, /„ oc 
v~ a , with additional widespread broad iron em ission and 
broad emission line s (IVanden Berk et all 1200 ll ). We fol- 
lowed iKellv et all (|2008l ) t o model the spectra, using 
an op tical Fe template from lVeron-Cettv. Jolv. fc Veron 
2004 ) and a UV Fe template from lVestergaard fc W ilkes 
2001). Each spectrum was simultaneously fit with a 



power-law continuum and an iron template using the 
Levenberg-Marquardt method for nonlinear mini- 
mization. We calculated the continuum luminosity from 
the power-law fit parameters. Light from the host galaxy 
can artificially inflate th e estimated AGN co ntinuum lu- 
minosity, as shown by iBentz et al.l (|2006h for the 35 
AGN with r everbe ration mapping data. In particular, 
IBentz et all ((2006) find that at 5100A, host galaxies 
contribute —20% of the measured flux for AGN with 
AL5100 > 10 44 erg/s, and -45% of the flux for AGN 
of A-L5100 < 10 44 erg/s. We note that the expected host 
contamination at L3000 or L1350 is much smaller, since 
most host galaxies (excepting very active star-forming 
hosts) have much less flux blueward of the 4000A break. 
So for the majority of our sample, the 150 AGN with 
masses measured from Mgn or Civ, we do not expect 
host contamination to be significant. However, the 15 
AGN with AL 300 o < 10 44 erg/s and H/3-derived masses 
may have luminosity estimates overestimated by a fac- 
tor of two, leading to black hole masses and Eddington 
ratios systematically overestimated by —0.15 dex. Fu- 
ture work in COSMOS will use host decompositions from 
HST/ACS images of z < 1 AGN to better characterize 
host contamination, but in this work we make no correc- 
tions for host galaxy flux. 

To fit the broad emission lines we subtracted the con- 
tinuum and Fe emission fits. Narrow absorption lines 
and narrow emission lines near the broad line (e.g., 
[Om]A4959 and [Om]A5007 near H/3) were fit by the 
sum of 1-2 Gaussian func tions and removed. Again fol- 
lowing EelBEjEai] (|2008h . each remaining broad emis- 
sion line profile was fit by the sum of 1-3 Gaussian 
functions, minimizing the Bayesian information crite- 
rion (BIC. ISchwartzl[l979h . Roughly, 2 Gaussians pro- 
vided the best fit for —40% of line profiles, while 1 or 
3 Gaussians each provided the best fit for —30% of line 
profiles. All fitting was interactive and inspected visu- 
ally, and if the multiple-Gaussian fit revealed a narrow 
(< 600 km/s) line in an H/3 emission line, the component 
was attributed to non-BLR origins and was removed. 
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Fig. 1. — Three spectra, representing fits and line profiles for 
Type 1 AGN with Civ (top), Mgll (middle), and H/3 (bottom). In 
the left panels, the spectra are shown by black histograms, and the 
dashed blue lines are the power-law plus iron emission continuum 
fits. At right, the blue dashed line shows the multiple-Gaussian fit 
to the continuum-subtracted line profiles, which are shown in black 
and shifted above for clarity. The fit to the H/3 line profile in the 
lower right panel includes a narrow (<r = 433 km/s) Gaussian which 
was removed before calculating the broad emission line width. The 
top two spectra were observed with Magellan/IMACS and have 
~10A resolution, while the bottom spectrum comes from the SDSS 
and has ~3A resolution. 

The FWHM was calculated directly from the multiple- 
Gaussian fit in order to minimize the effects of noise in 
the original spectra. The multiple-Gaussian fit was ro- 
bust to a variety of line profiles, and simulated spec- 
tra (see §3.2) revealed measured FWHM errors of only 
&FWHM /FWHM ~ 10%. 

Three examples of spectra with fitted continua and 
continuum-subtracted line profiles are shown in Figure 
[TJ The spectra are representative of typical fits for 
each of the Civ, Mgll, and H/3 emission lines. At left 
the power-law and iron emission fits are shown by the 
dashed blue lines. The right panel shows the multiple- 
Gaussian line profile fits as dashed blue lines, with the 
continuum-subtracted line profile shifted above by an ar- 
bitrary amount for clarity. The fit to the H/3 line profile 
in the bottom right panel includes a narrow (a = 433 
km/s) Gaussian which is not associated with the BLR 
and was removed. Even for noisy spectra like the middle 
panel, the spectral fitting provides a robust continuum 
and isolates the emission line. 

3. ESTIMATED BLACK HOLE MASSES 

We estimate black hole masses using our mea- 
sured bro ad line velocity w i dths and the scaling re- 
lations of IVestergaard et alj (|in prep.1 ) for Mgll and 
IVestergaard fc Peterson! ([2006T ) for H/3 and Civ. These 
relations all take the form of Equation 1, with XL\ in 
units of 10 44 erg/s and Vfwhm hi units of 1000 km/s; 
A = 6.91, B = 0.50, and A = 5100A for H/3; A = 6.86, 
B = 0.50, and A = 3000A for Mgll; A = 6.66, B = 0.53, 
and A = 1350A for Civ. The Mgll relation was derived 
from SDSS quasars with both Civ and Mgll in the spec- 
trum, and it is designed to produce black hole masses 
consistent with those measured from Civ. In our sam- 
ple, we measure H/3 for 32 AGN, Mgn for 134 AGN, and 
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Fig. 2. — Black hole mass versus optical luminosity and redshift 
for 182 Type 1 AGN in COSMOS. The black hole masses are de- 
rived from the broad line velocity widths of H/3 (blue diamonds), 
Mgll (green di amonds), and Civ (red diamonds) depending 
on the redshift (Vestcreaard & Peterson 2006; Vcstcrgaard ct al. 
I in prep]) . The black hole masses have uncertainties o f ~0.4 dex 
from the scaling relations (Krolik 2001: Shcn ct al. 2008), as shown 
by the typical error bars are shown in the lower right of each 
panel (the redshift error is < 1%). The diagonal tracks at left 
represent Eddington ratios, assuming a bolometric correction of 
Lbol = 5(AL3ooo) (Richards ct al. 2006J). At rig ht the gray crosse s 
show scaling relation masses of SDSS quasars (|Kellv et al.ll2008h . 
The depth of COSMOS allows us to probe low mass and weakly 
accreting SMBHs, revealing that L/L Edd < 0.01 Type 1 AGN do 
not exist. 

Civ for 38 AGN (19 SDSS AGN have both Mgn and Civ, 
and 3 have both H/3 and Mgn). AGN with estimates of 
Mbh from two different emission lines are treated as two 
separate objects in our subsequent analyses. 

Table Q] presents the catalog of black hole masses and 
line measurements. AGN with both Mgn and Civ or 
both H/3 and Mgn present have two entries in Table [I] 
one for each emission line. The full catalog then contains 
204 entries: 182 Type 1 AGN with M BH estimates, 22 of 
which have two sets of broad emission line measurements. 
The black hole masses are shown with continuum lumi- 
nosity (calculated from the power-law fit) and redshift 
in Figure O The diagonal tracks in the figure represent 
Eddingt on ratios using a bo lometric correction of 5 for 
AL3000 ([Richards et alj [20061 ) . We als o show a compari- 
son sample of brighter SDSS quasars (jKellv et al.ll2008f l 
in order to highlight the lower black hole masses probed 
by COSMOS. 

3.1. Error 

The scaling relations have uncertainties of ^0.4 dex, 
altho ugh there may be larger syst e matic uncertain- 
ties (iKrolikl [27)01 [Collin et al.1 12001 iFine et all [20081 : 
iMarconi et al.ll200g ). Measurement errors in the lumi- 
nosity and emission line FWHM also contribute, but the 
uncertainty from the scaling relations dominates. We 
test the luminosity error in the left panel of Figure [31 
which compares the luminosit y estimates from this wor k 
to duplicate estimates frorn_JVierloni et al.l (jin nreri.fi . 
The luminosity estimates of M erloni et al.l (jin prep.1 ) use 
independent redshifts from VLT/VIMOS spectra and are 
calculated from a fit to the IR to X-ray multiwavelength 
spectral energy distribution, instead of from the opti- 
cal spectrum itself (as in this work). The scatter be- 
tween the two luminosity estimates is a = 0.25 dex. The 



4 



Trump et al. 



TABLE 1 

COSMOS Type 1 AGN Black Hole Mass Catalog 



Object 


Rcdshift 


S/N 


loe(XL • ) 

by 3000A' 


iu^ajq .5 — 2kcV ) 


Line 


FWHM 




(J2000) 




(per pixel) a 


[erg/s] 


[erg/s] b 




(km/s) 




SDSS J095728. 34+022542.2 


1.54 


7.0 


45.00 


44.47 


Mgll 


4491 


8.665 


SDSS J095728. 34+022542. 2 


1.54 


7.0 


45.00 


44.47 


CIV 


2776 


8.135 


COSMOS J095740.78+020207.9 


1.48 


17.9 


43.46 


44.41 


Mgll 


6701 


8.244 


SDSS J095743. 33+024823.8 


1.36 


3.4 


44.60 


43.50 


Mgll 


3472 


8.243 


COSMOS J095752. 17+015120.1 


4.17 


7.3 


45.19 


44.36 


CIV 


4603 


8.656 


COSMOS J095753. 49+024736.1 


3.61 


4.8 


44.14 


44.53 


CIV 


2629 


7.997 


SDSS J095754. 11+025508.4 


1.57 


6.1 


45.07 


44.45 


Mgll 


4500 


8.701 


SDSS J095754.70+023832.9 


1.60 


8.0 


45.17 


43.69 


Mgll 


3361 


8.498 


SDSS J095754.70+023832.9 


1.60 


8.0 


45.17 


43.69 


CIV 


6384 


8.946 


SDSS J095755. 08+024806.6 


1.11 


8.7 


44.92 


44.08 


Mgll 


3574 


8.426 


COSMOS J095755. 34+022510.9 


2.74 


3.2 


44.38 


-1.00 


CIV 


3879 


8.074 


COSMOS J095755.48+022401.1 


3.10 


19.9 


45.28 


45.00 


CIV 


3527 


8.436 



a The SDSS spectra have 3 pixels per resolution element, and the Magellan/IMACS spectra have 5 pixels per resolution 
element. b AGN with no soft X-ray detection have an entry of -1.00 for log(Lo.5-2kev)- 



Type 1 AGN in COSMOS have an average variability 
of -0.15 dex (|Salvato et al.ll2009t) . the remaining lumi- 
nosity scatter can be attributed to the different methods 
of estimates. Since Mbh ~ £ ' 5 , our luminosity error 
contributes very little to the overall Mbh uncertainty. 

Line measurements of synthetic spectra (described 
in §3.2 below) show that our FWHM error is only 
ctfwhm/FWHM ~ 10% at f£ B - 22. The right panel 
of Figure [3] compares the duplicate estimates of Mbh for 
spectra with two broad emission lines. Red diamonds in- 
dicate spectra with both Mgll and H/3, while blue crosses 
indicate both Mgll and Civ. The scatter between the 
different estimates is only a = 0.36 dex, nearly the same 
as the expected intrinsic scatter for Mbh- This sug- 
gests that the statistical error in the mass estimators is 
not correlated to the choice of emission line (see also 
iKellv fc Bechto ld 2007). If there were systematic offsets 
in the mass estimators, they would cause a constant shift 
in the mass estimate for each line, and therefore would 
not contribute to the scatter between two lines. The sta- 
tistical intrinsic scatter, however, would not "cancel" in 
such a way. Because the scatter between lines is compa- 
rable to the expected intrinsic scatter, our estimates of 
Mbh probably do not have significant systematic errors. 

3.2. Completeness 

.Previous work (jTrump et alj 12009ft tested the com- 
pleteness of the Type 1 AGN sample, with simulated 
spectra showing that Type 1 AGN are correctly iden- 
tified (with high confidence redshifts) at 90% complete- 
ness to S/N>2.87. But even for correctly identified Type 
1 AGN, our measurements of Mbh are roughly limited 
by spectral S/N and FWHM, since we cannot identify 
or measure broad emission lines for spectra that have 
lines so broad that they become confused with noise or 
the Fe emission. To test these limits, we create 100 
synthetic spectra of Type 1 AGN with Mgll in the ob- 
served wavelength range. We choose Mgll because it 
is the most common line used for our Mbh calcula- 
tions and also because it is the broad emission line most 
contaminated by widespread iron emission. These syn- 
thetic spectra are formed by making a composite of all 
1 < z < 2.4 observed spectra, then removing the rest- 
frame 27001 < A < 2900A Mgll region. To each of the 
100 spectra we then re-add a Mgll region with random 



44.5 



s 44.0 



43.5 
1.0 
0.5 
0.0 
-0.5 
-1.0 



9.5 



43.5 44.0 44.5 45.0 45.5 
VIMOS log(XL 3000 ) [erg/s] 



7.5 
1 .0 
0.5 
0.0 
-0.5 
-1.0 
7 



■*■+ +/ 

V + 



o o"/ + 



• , T 

+ + + + + 



5 8.0 8.5 9.0 9.5 10.0 
Hb or CIV log(M BH ) [M„„] 



Fig. 3. — At left, the luminosity estimates for 48 AGN in 
this work are compare d to estimates of the same objects from 
IMerloni et al.l din prep.1) , which have independent redshifts from 
VLT/VIMOS spectra and are calculated from a fit to the multi- 
wavelength SED. The scatter between the two luminosity estimates 
is <t = 0.25 dex, which contributes very little to the overall scatter 
in Mbh- At right, Mbh estimates are compared for the 22 spec- 
tra with two emission lines (either Mgll and H/3, red diamonds, or 
Mgll and Civ, blue crosses). The scatter between separate Mbh 
estimates is only a = 0.36 dex, indicative that the intrinsic scatter 
of ~0.4 dex dominates the Mbh error. 

FWHMs and line areas, and then each spectrum has ran- 
dom noise added. We choose the line areas and noise to 
be normally distributed in the ranges of measured line 
area and S /N in the original spectra, while the FWHMs 
are chosen to probe our sensitivity to the broadest line 
widths, FWHM> 10000 km/s. 

We show the FWHM and S/N of both the observed and 
simulated spectra in Figure 2] The blue diamonds and 
red squares represent the simulated spectra that we suc- 
cessfully measure and those we miss, respectively. The 
vertical dotted line shows the sample's 90% complete- 
ness limit for correctly identifying high-confidence Type 
1 AGN. The dashed line shows the limit to where we can 
measure broad emission lines, corresponding to > 95% 
completeness since only one "missed" synthetic spectrum 
lies below the line. The observed Type 1 AGN tail off well 
before the dashed line. Translating FWHM and mag- 
nitude into black hole mass, as an example, a Type 1 



AGN with i 
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22 at z — 1.7 with a Mgll profile of 
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Fig. 4.— The broad emission line FWHM and S/N for our 182 
Type 1 AGN are shown as black crosses. The upper axis shows the 
*AB magnitude scale which roughly corresponds to the S/N (from 
Figure 11 of ITrump eTaTI |2009|) . Also shown are the measured 
FWHM and S/N for 100 simulated Type 1 AGN spectra: blue di- 
amonds represent spectra we can correctly identify and measure, 
while red squares represent synthetic spectra that we incorrectly 
identify or cannot properly measure FWHM. The dashed line shows 
the limit to which we correctly identify and measure > 95% of syn- 
thetic spectra. The dotted vertical line shows the 90% complete- 
ness to identifying and assigning high- confidence redshif ts to Type 
1 AGN, from additional simulations (Trump ct al. 2009). The lack 
of high-FWHM AGN identified in COSMOS is not from selection 
effects, since we can correctly identify and measure synthetic Type 
1 AGN with much broader lines at S/N> 3. 

FWHM= 10000 km/s would have log(Af B #/M ) - 10.6 
and L/LEdd ~ 0.0002. We do not detect such objects in 
our sample, yet our simulations show that they do not 
lie beyond our detection limits. 

4. DISCUSSION 

In Figure [2] all Type 1 AGN lie within the re- 
gion of 0.01 < L/LE dd < 1, a result supported by 
iKollmeier et al.l ((2006) . This implies that the broad 
emission line region of Type 1 AGN might become un- 
detectable as the accretion drops below L ~ O.Olisdd- 
Such objects might be observed as unobscured Type 2 
AGN, the possible remnants of "dead" Type 1 AGN 
whose accretion dis k geometries changed as their accre- 
tion rates fell (e.g.. iHopkins et ail [2008). Or these low 
accretion rate AGN may be diluted, with their emission 
falling below the light of their host galaxy. The AGN 
emission may also be unable to blow out local obscur- 
ing material, causing their BLR to lie undetected behind 
obscuration. 

To study the accretion rates of our AGN, we calculate 
the intrinsic luminosity from o ur measured ALqqoo a nd 
£o.5-2keV, using the relations of iMarconi et alj ((2004): 

log[L 7 /(A£ 3 ooo)] = 0.65 - 0.067£ + 0.017£ 2 - 0.0023£ 3 

(2) 

log(L//L . 5 _ 2 ) = 1.65 + 0.22£ + 0.012£ 2 -0.0015£ 3 (3) 

Here C = log(Li) —45.58 and all luminosities are in units 
of erg/s. The intrinsic luminosity Lj is designed to be 
a bolometric luminosity which excludes reprocessed (IR) 
emission, so that the Eddington ratio Lj / L Edd represents 
a robust measure of the accretion onto the black hole. 
We use the Newton method to solve each equation for 
Lj from AL3000 and £0.5-2- We then average the two 



values of Lj for our final value (excepting 7 AGN where 
we estimate Li from AL3000 only because they lack soft 
X-ray detections). 

We show the Eddington ratios of our Type 1 AGN 
with Mbh 7 AL3000 ) and redshift in Figure [5l The dia- 
monds show individual objects, while the solid lines show 
the means and scatter in equal-sized bins. The scatter 
(the standard deviation of the mean) is calculated as the 
square root of the second moment of the data in each 
bin. Our mean Eddington ratio for all Type 1 AGN 
is Lj/LEdd ~ 0.1, lower than th e value of LjlL^dd ~ 
0.3 found in previous surveys (jKollmeier et al.l 120061 : 
iGavignaud et alJl2008ft. This is partly expl ained by the 
depth of COSMOS~ lKollmeier et aD (f2006h noted that 
their R < 21.5 AGES sample was only complete to 
LijLEdd ~ 0.1, while our simulations in §3.2 show that 
COSMOS can reach much weaker accretors. In addition, 
the center panel of Figure [5] shows that accretion rate 
may increase with optical/UV luminosity, suggesting an 
additional reason for our lower m ean Eddington ratio: 
most o f the IKollmeier et all (|2006f) and IGavignaud et al.1 
((2001 AGN have AL3000 > 10 45 erg/s, where our AGN 
have LijLEdd ~ 0.2. The majority of our AGN have 
A£3ooo < 10 45 erg/s and so we find a lower mean accre- 
tion rate. 

The apparent decrease in accretion rate with black hole 
mass and the apparent increase in accretion rate with 
redshift can be explained by selection effects: low accre- 
tion rate AGN are more difficult to detect if they are also 
low mass or at higher redshift. The increase in accretion 
rate with optical l umino sity, however, is also observed by 
IGavignaud et al. (2008) and is probably a physical effect. 
We performed a linear regression analysis of the corre- 
lation between accretion rate and optical luminosity us- 
ing the publ icly available IDL program linmix_err .pro 
l|Kefi^l2007h . Using errors of 0.25 dex in log(AL3ooo) 
and 0.4 dex in log(L/ / LEdd)-, linear regression indicates 
that log{L I /L Ed d) ~ (0.28 ± 0.06) log(AL 3000 ). In other 
words, accretion rate is correlated with opticaly lumi- 
nosity at the 4.8<7 level. As a Type 1 AGN increases 
in accretion rate, its optical emission becomes a larger 
fraction of its total bolometric output because its cool 
accretion disk emi ts more brightly. This is consistent 
with the results of Ke lly et "alT(|2008t ). which show that 
®ox (the ratio between optical/UV and X-ray flux) be- 
comes more X-ray quiet with accretion rate. Thus a more 
rapidly accreting Type 1 AGN has more of its emission 
in its cool (optical) disk than in its hot (X-ray) corona, 
possibly because the disk grows larger or thicker as the 
accretion rate approaches the Eddington limit. 

The scatter (square root of the second moment) in 
Li/LEdd, shown in the bottom panels of Figure is 
typically only ~0.4 dex in each bi n. This is greater than 
previously measured d i spersions (IKollmeier et alj 120061 : 
IGavignaud et al.l 120081 : iFine et aLll2008f) . and indicates 
that COSMOS is more sensitive to low accretion rate 
Type 1 AGN than previous studies. Yet it is remarkable 
that the dispersion is not larger than the scatter from 
the scaling relations: the intrinsic dispersion in Edding- 
ton ratio might then be ~0, with nearly all Type 1 AGN 
of a given mass, luminosity, and/or redshift a ccreting at 
a very narrow range of of accretion rates. IFine et aLl 
(2008) note that at such low measured dispersions, if the 
intrinsic dispersion in accretion rate is much greater than 
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Fig. 5. — The Eddington ratio (accretion rate) with black hole mass, luminosity, and redshift for our Type 1 AGN. Eddingt on ratio was 
calculated using an intrinsic luminosity Lj estimated from AL3000 and £o.5-2kcV using the relations of Marconi ct al. (2004). Diamonds 
represent individual objects with masses estimated from H/3 (blue), Mgll (green), and Civ (red). The large crosses in the top plots show 
the mean accretion rate in each bin of Mbh or redshift, while the gray lines at the bottom show the standard deviation (the square root 
of the second moment) in each bin. The dispersion deviation is also shown by the vertical error bar in the top plots. Bins were chosen to 
each have the same number of objects. Selection effects cause the apparent trends of Lj /L^dd decreasing with Mbh and increasing with 
redshift, but the increase of Lj /L^dd with AL3000 is a physical effect caused by changes in the accretion disk. The dispersion is generally 
~0.4 dex, higher than that of previous, less sensitive surveys. 



0, then the BLR cannot be in a simple virial orbit. Ac- 
curate scaling relations would th en require a luminos ity- 
dependent ionization parameter ( Marconi et al.ll2008f) or 
a more complex BLR geometry ()Fine et al.l 12008). We 
note, however, that the myriad uncertainties involved 
in estimating Lj/LEdd make concrete conlusions diffi- 
cult. And although we find significant evidence that 
Lj/LEdd < 0.01 Type 1 AGN do not exist (or are very 
rare), the spectroscopic flux limit may still miss some 
Lj/LEdd ~ 0.01 AGN at lower luminosities. 

5. SUMMARY 

The black hole masses of Type 1 AGN in COSMOS in- 
dicate that Type 1 AGN accrete at a narrow range of high 
efficiencies, Lj/Lsdd ~ 0.1. When the accretion rate of 
an AGN lowers, less of its luminosity is emitted opti- 
cally. When a Type 1 AGN accretion rate drops below 
Li/LEdd ~ 0.01 the BLR becomes invisible, due to ob- 
scuration, dilution, or an altered accretion disk geometry. 
We additionally measure higher dispersions in accretion 
rate than previous, less sensitive surveys, although the 
dispersion is still no lar ger than the i ntrins ic uncertainty 
in the scaling relations iKellv et al.l (|2008f ) find that the 



bolometric correction depen ds on black hole mass, and 
I Vasudevan fc F abian (2009) find that it correlates with 
Eddington ratio. This makes characterizing the distri- 
butions of Li/LEdd and its scatter rather difficult. We 
partially mitigate the systematic uncertainties by using 
both AL3000 and Lo.5-2 to estimate Lj. Future work in 
COSMOS will use more accurate bolometric luminosities 
calculated from the full multiwavelength dataset. 



We thank Brad Peterson for helpful discussions on 
emission line measurements and Misty Bentz for com- 
ments on host contamination. We thank Marianne 
Vestergaard for providing the Fe templates and for dis- 
cussions on the scaling relations. Amy Stutz and Aleks 
Diamond-Stanic provided useful discussions on statis- 
tics. JRT acknowledges support from NSF ADP grant 
NNX08AJ28G and an ARCS fellowship. BCK acknowl- 
edges support from NASA through Hubble Fellowship 
grant #HF-01220.01 awarded by the Space Telescope 
Science Institute, which is operated by the Association of 
Universities for Research in Astronomy, Inc., for NASA, 
under contract NAS 5-26555. 



REFERENCES 



Bentz, M. C, Peterson, B. M., Pogge, R. W., Vestergaard, M., & 

Onken, 2006, C. A. ApJ, 644, 133 
Bentz, M. C, Peterson, B. M., Pogge, R. W., & Vestergaard, M. 

2008, ApJ, 694, 166 
Bigelow, B. C, Dressier, A. M., Shectman, S. A., & Epps, H. W. 

1998, in Proc. SPIE Vol. 3355, 225, Optical Astronomical 

Instrumentation, Sandro D'Odorico; Ed. 
Bongiorno, A. et al. 2007, A&A, 472, 443. 
Brandt, W. N., & Hasinger, G. 2005, ARA&A, 43, 827. 
Brusa, M. et al. 2007, ApJS, 172, 353 
Brusa, M. et al. 2009, ApJ, in preparation 
Cappelluti, N. et al. 2008, A&A submitted 

Collin, S., Kawaguchi, T., Peterson, B. M. & Vestergaard, M. 2006, 
A&A, 456, 75 



Davies, R. I. et al. 2006, ApJ, 646, 754 

Di Matteo, T., Springel, V., & Hernquist, L. 2005, Nature, 433, 
604 

Dietrich, M. & Hamann, F. 2004, ApJ, 611, 761 
Ferrarese, L. & Merritt, D. 2000, ApJ, 593, 9 
Fine, S. et al. 2008, MNRAS, 390, 1413 
Gavignaud, I. et al. 2008, A&A, 492, 637 
Gebhardt, K. et al. 2000, ApJ, 539, 13 
Greene, J. E., & Ho, L. C. 2006, ApJ, 641, 21 

Hopkins, P. F., Hernquist, L., Cox, T. J., Di Matteo, T., Robertson, 

B., & Springel, V. 2006, ApJS, 163, 1 
Hopkins, P. F., & Hernquist, L. 2006, ApJS, 166, 1 
Hopkins, P. F., Hickox, R., Quataert, E., & Hernquist, L. 2008, 

MNRAS submitted (arXiv/0901.2936) 



COSMOS Black Hole Masses 



7 



Kaspi, S., Smith, P. S., Netzer, H., Maoz, D., Jannuzi, B. T. & 

Giveon, U. 2000, ApJ, 533, 631 
Kaspi, S., Brandt, W. N., Maoz, D., Netzer, H., Schneider, D. P. 

& Shemmer, O. 2007, ApJ, 659, 997 
Kelly, B.C. 2007, ApJ, 665, 1489 
Kelly, B. C. & Bechtold, J. 2007, ApJS, 168, 1 

Kelly, B. C, Bechtold, J., Trump, J. R., Vestergaard, M., & 

Siemiginowdka, A. 2008, ApJS, 176, 3557 
Koekemoer, A. M. et al. 2007, ApJS, 172, 196 
Kollmeier, J. A. et al. 2006, ApJ, 648, 128 
Kormendy, J. & Richstone, D. 1995, ARA&A, 33, 581 
Krolik, J. H. 2001, ApJ, 551, 72 
Magorrian, J. et al. 1998, AJ, 115, 2285 

Marconi, A., Risaliti, G., Gilli, R., Hunt, L. K., Maiolino, R. & 

Salvati, M. 2004, MNRAS, 351, 169 
Marconi, A. et al. 2008, ApJ, 678, 693 
McLure, R. J., & Jarvis, M. J. 2002, MNRAS, 337, 109 
Merloni, A., Bongiorno, A., Trump, J. R. et al. in prep. 
Onken, C. A., Fcrrarese, L., Merritt, D., Peterson, B. M., Pogge, 

R. W., Vestergaard, M., & Wandel, A. 2004, ApJ, 615, 645 
Onken, C. A. et al. 2007, ApJ, 670, 105 
Peterson, B. M. & Bentz, M. C. 2006, NewAR, 50, 796 
Richstone, D. et al. 1998, Nature, 395, 14 



Richards, G. T. et al. 2006, ApJ, 166, 470 
Salvato, M. et al.2009, ApJ, 690, 1250 
Schneider, D. P. et al. 2007, AJ, 130, 367 
Schwartz, G. 1979, Ann. Statist., 6, 461 
Scoville, N. et al. 2007, ApJS, 172, 38 

Shen, Y., Greene, J. E., Strauss, M. A., Richards, G. T., & 

Schneider, D. P. 2008, ApJ, 680, 169 
Silk, J. & Rees, M. J. 1998, A&A, 331, 1 
Soltan, A. 1982, MNRAS, 200, 115 
Trump, J. R. et al. 2009, ApJ, 696, 1195 
Ueda, Y. et al. 2003, ApJ, 598, 886 
Vanden Berk, D. E. et al. 2001, AJ, 122, 549 
Vasudevan, R. V. & Fabian, A. C. 2009, MNRAS, 392, 1124 
Veron-Cetty, M.-P., Joly, M., & Veron, P. 2004, A&A, 417, 515 
Vestergaard, M. & Wilkes, B. J. 2001, ApJS, 134, 1 
Vestergaard, M. 2004, ApJ, 601, 676 
Vestergaard, M. & Peterson, B. M. 2006, ApJ, 641, 689 
Vestergaard, M. et al. in prep. 

Wu, X.-B., Wang, R., Kong, M. Z., Liu, F. K., & Han, J. L. 2004, 

A&A, 424, 793 
York, D. G. et al. 2000, AJ, 120, 1579 



